Diverging magnetothermal response in the one-dimensional Heisenberg chain 
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A current of magnetic moments will flow in the spin-1/2 Heisenberg chain in the presence of 
an external magnetic field B and a temperature gradient AT along the chain. We show that this 
magnetothermal effect is strictly infinite for the integrable Heisenberg-model in one dimension. 
We set-up the response formalism and derive several new generalized Einstein relations for this 
magnetothermal effect which vanishes in the absence of an external magnetic field. We estimate 
the size of the magnetothermal response by exact diagonalization and Quantum Monte Carlo and 
make contact with recent transport measurements for the one- dimensional Heisenberg compound 
Sr 2 Cu0 3 . 

PACS numbers: 75.30.Gw,75.10. Jm,78.30.-j 



I. INTRODUCTION 

The nature of magnetic and thermal transport in mag- 
netic insulators with reduced dimensions is a long stand- 
ing problem. Huber, in one of the first works on the 
subject j* evaluated the thermal conductivity n(T) for 
the Heisenberg-chain with an equation-of-motion approx- 
imation and found a finite n(T), a result which is, by 
now, known to be wrong. A few years later, Niemeijer 
and van Vianen calculated n(T) in the case J z — us- 
ing the Jordan- Wigner transform and found a diverging 
result £ It has been shown recently^ that the energy- 
current operator commutes with the Hamiltonian for the 
spin-1/2 Heisenberg chain. The thermal conductivity 
is consequently infinite for this model. The intriguing 
question, "under which circumstances does an interacting 
quantum-system show an infinite thermal conductivity" , 
is being intensively studied theoretically^^ motivated, 
in part, by new experimental findings. 

An anomalous large magnetic contribution to 
k has been observed for the spin-ladder system 
Sri4_ x Ca x Cu2404i and CagLa5Cu2404iji& raising the 
possibility of ballistic magnetic transport limited only 
by residual spin-phonon and impurity scattering. Large 
energy-relaxation times have been found in recent 
experimental and theoretical studies of the thermal 
conductivity for the quasi-one dimensional spin-chain 
compounds SrCuC>2 and S^CuOs, consistent with 17 O 
NMR studiesJfi 

It is knownj&ii that there is no magnetothermal ef- 
fect in the Heisenberg chain since the magnetic-current 
operator is a pseudo- vector. This can be seen by a sim- 
ple symmetry argument in the standard setting, when we 
consider only a temperature gradient along the sample. 
A non-zero magnetothermal effect would yield a magneti- 
zation along an arbitrary quantization axis. But because 
of the isotropic conditions there is no preferred direction 
along which such a magnetization might occur, the effect 
vanishes consequently. 

The situation is, however, different if there is an ex- 
ternal magnetic fieldi A temperature gradient will now 
cause a magnetization current with a magnetization vec- 



tor parallel to the field, see Fig. ^ This magnetother- 
mal effect can be considered the generalization to mag- 
netic systems of the Secbcck effect occurring in nor- 
mal conductors. It is diverging for the isolated Heisen- 
berg chain, leading to a finite non-zero 'magnetothermal 
Drude weight'. In this paper we discuss magnetothermal 
effects in spin chains, i.e., we examine how the magnetic 
and thermal currents couple to external sources. 

In Sec. [n] we will discuss the differences in between 
the thermomagnetic effects found in normal metals and 
the magnetothermal response we study here for magnetic 
insulators. In Sec. IIIII IIVI and in Sec. we set-up the 
formalism for the magnetic and energy current operators 
in the presence of an external magnetic field and for the 
correlation and response- functions, respectively. 

In Sec. IVII we evaluate the magnetothermal response 
for the xy-chain via a Jordan- Wigner transformation and 
derive and discuss several generalized Einstein relations 
in Sec. IVHI In Sec. IVIIII we derive and discuss a new 
exact magnetothermal Einstein relation. 

We present in Sec. IIXI numerical results obtained 
from Quantum-Monte-Carlo Simulations and exact- 
diagonalization studies. Finally, we discuss the size of the 
magnetothermal effect expected for S^CuOa together 
with a dimensional analysis in Sec. Ixl 



II. THERMOMAGNETIC AND 
MAGNETOTHERMAL RESPONSE 



In a conventional solid there is normally a variety of 
ways in which external parameters (e.g., current density 
j, temperature gradient and magnetic field B) couple to 
the voltage V/x and the thermal current je leading to 
numerous thermoelectric and thermomagnetic effects. A 
phenomenological equation which gives full account of all 
these effects would look as (see, e.g., Ref. IT^) 
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(Bx is understood to be applied to both components 
separately.) p is the resistivity, k the thermal conduc- 
tivity, and the other coefficients describe the Seebeck- 
(Q), Peltier- (II), Hall- (R), Nernst- (TV), Ettingshausen- 
(E), and Righi-Leduc-effect (L). The latter four effects 
are called thermomagnetic - for obvious reason - and 
are not to be confounded with the effects which are to 
be discussed in this paper which we denote magnetother- 
mal effects. One should note that in magnetic systems 
the magnetic field corresponds to the chemical potential, 
leaving no analo gues for the entries of the second matrix. 

Following Ref. 12|, we write the dependence of the cur- 
rent densities jm (particle) and je (energy) on external 
sources (field VB = — V/i and temperature gradient VT) 
in the following manner 



Jm 
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Lmm 
Lem 



Lme 
Lee 



VB \ 
(-VT) ) 



(1) 



[We deviate from the standard notation by indexing 
the components of the tensor by M and E rather than 1 
and 2. The notation with M as particle index and B as 
chemical potential borrows from spin systems, this arti- 
cle's chief case of interest. Eq. may therefore be used 
as a reference for both electrical and magnetic systems.] 
Normally, the two off-diagonal components governing the 
Seebeck and Peltier effect are not independent but related 
via the Onsager relation. Here, with the given choice of 
external forces, it reads (where T denotes temperature) 
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To compute the response functions L\j we use the stan- 
dard Kubo-formula2iSiiii4 which yields for the Heisen- 
berg model generally an L- Tensor with infinite compo- 
nents, i.e., a delta function times a weight- factor. These 
weight-factors will be denoted by the entries of the L- 
Tensor Lij (without "-accent). Normally, one has an 
additional finite contribution - the regular part - which 
will turn out to be zero for all but the MM-component. 
We assume that due to some - so far unaccounted - scat- 
tering processes with external degrees of freedoms (like 
phonons or impurities) the infinite peak broadens, and 
the coefficients (we consider only the real part) may be 
replaced according to (cf., e.g., Ref. 14) 



Lij(u>) = Lij7r5(u>) 



LijT 



1 + {uJTf 



(3) 



with i,j G {M,E} and a finite relaxation time r - for 
simplicity, a possible dependence on i and j shall be ne- 
glected. 



standard xirz-chain, H = VJ n h n with 
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Hb is the Bohr magneton and the lattice spacing will be 
denoted by c ; the S n are supposed to be dimensionless. 

Instead of the electrical current we consider the mag- 
netic current 



c J 



(5) 



= [S„ x S„ + i] 2 = i — [H, Pm] > 

n 

where we have defined the magnetic polarization Pm = 
Hb Yin n ^rn which is possible only for chains with open 
boundary conditions (OBC). 

For the energy currenl££ we have equivalently 

— Je = i ^ [h n , h n+ i] — B — Jm + ^ S ra ■ (S ra+ i x S n+2 ) 

n n 

= B-J M +J2 dct [S„S„ +1 S„ +2 1 =i[H,P E ], (6) 



where we have used S„ = (S%, S^, S„) and the definition 
S,i = {J x Sn, JySv, J Z S Z ). In the following we will denote 
with ji = J; /Vol, i € {M,E} the respective current 
densities. Here Vol = cN is the one-dimensional volume, 
where c is the lattice constant and V the number of sites. 
The energy polarization Pe = Yn entering Eq. 
for OBC represents a temperature gradient, and, 
as is well known^& entails always a magnetic polariza- 
tion Pm, which stems from the chemical- potential-term 
(magnetic- field-term). To see this we insert a site- 
dependent (linear) 

13 = (3{n) =f3 + ncVP 

into the thermodynamical expectation valued The 
Boltzmann- factor can be rewritten as follows: 

exp ^-^2f3(n)h^j = exp -/3 + nh n 



cxp [-0(H + FcP E )] 



(7) 



where F = V/3//3 = -(VT)/T. This motivates the use 
of Pe to model a temperature gradient. 



IV. CORRELATION FUNCTIONS 



III. CURRENT OPERATORS 



Setting up the notation for the response theory we de- 
fine with 13 ; 17 



In this article, we discuss the analogue of the ther- 
moelectric effect in spin-chains. The Hamiltonian is the 



A(AB)(z) ee - 



e izt ([A(t),B])dt (8) 
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the retarded Green's function of two operators A and 
B. In our case A and B will be mostly current opera- 
tors; we therefore introduce the notation A(JjJj) = A{j. 
A(AB)(z) vanishes when one of the operators is a con- 
stant of motion and commutes with the Boltzmann- 
factor: by the cyclic property of the trace one has in 
this case {AB) = (BA). 

The energy-current operator Je commutes with the 
Hamiltonian in the absence of an external magnetic field 
B. Using Eq. J6J we have A ME — A E m = BA MM and 
A-ee = B 2 A MM . 

The isothermal susceptibility!^ is given by 



X T (AB) = f(AA(r)AB)dr , 
Jo 



(9) 



with AA = A- (A) and AB = B-(B). It takes the usual 
form x T (AB) = (3(AAAB) when one of the operators 
is a constant of motion. The generalized Drude weight 
(which is used, e.g., in Ref. 4J) is defined by 



/>oo 

((AB)) = lim(-iz) / e ltz (AA(t)AB) dt 



(10) 



These three correlation functions are not 
independent, comparing the respective 

eigenstates-representation, e.g., ZA(AB)(0) = 
J2e M ^A nk AB kn (e-P E * - e-^«) / (E n - E k ), 
where Z is the partition function, one finds 



f3 ((AB)) + A(AB)(0) = X T (AB) 



ill) 



Under open boundary conditions, when the relations 



[H, Pj] are valid, we can express the isothermal 



susceptibility, Eq. 0, 



X 



T ' OBC {J^) 



(12) 



as a static expectation value. Note, however, that the 
value of the susceptibility is independent of the boundary 
conditions in the thermodynamic limit: X T (JiJj) — 



V. MAGNETOTHERMAL COEFFICIENTS 

We now discuss the general recipe for the computation 
of the entries of the L-Tensor, appearing in Eq. ©. 

We assume a perturbation in the form of a polariza- 
tion, (cf., e.g., Ref. 13) i.e., we add to the Hamiltonian 
a term cVB ■ Pm with Pm '■= I^b J2 n n ^n an< ^ compute 
the response of the current density operator to obtain the 
MM-component. (A different but equivalent approach 
is given in Ref. Il7l) To compute the remaining entries 
of the L-Tensor, we replace the magnetic by the energy 
current density (in the second row) and/or substitute 
Pe ■= J2 n n ^ n f° r P M ~ as well as cVT for cVB - and 
add a factor ks{3 (in the second column). [One should 
compare this procedure with Eq. Q).] 



The linear response theory gives the following contri- 
bution to Lij-i& 



(k B py~ l l / e^ifaPj (~t)])dt 



where (3 M := (3, (3 E := /3 2 , etc., and j. t = J 4 /Vol (see Sec. 
IIII|) . Integration by parts yields for the above expression 



lim^- 

z^o — iz 



" v ' 



One should note that Ajj(O) = if one of the currents 
commutes with the Hamiltonian, and consequently there 
is no regular part. The factor (— iz)^ 1 accounts for the 
delta-function in the zero-frequency limit and we have 



L; 



It is also possible to express this result by the correlation 
function defined in Eq. ljTU|) using Eqs. ljTT|) and lfT^|) : 



La = (kepy- 1 [x T (jiJj)-kj(P)] = (&■/,■» ■ 

; (14) 

From Eq. (|13|) follows the well-known result for the mag- 
netic Drude weight Lmm = (—7mm) § — Amm, where 
Tmm — -'[jm, Pm] is the kinetic energy per site - apart 
from a prefactor. In view of Eq. (|14J) we define the gen- 
eralized kinetic energy 



T 



-i\j i (B = 0),P j (B = 0)}. 



(15) 



Here, one should note a particularity produced by the 
twofold role played by the operator Pe in the xxz-chain: 
it not only describes thermal response but also acts as 
a 'boost '-operator for the constants of motion^ Hence, 
Je as well as Tee are constants of motion. 

Taking into account that Je commutes with the Hamil- 
tonian for B — and using Eqs. (jBJ) and (|13fl we have 

Lem = «t([je, Pm]) - BA M m 
h 

= i^(\j E (B = 0),P M ]) + B{i^( [ JM , Pm ] ) - A M m } 

= (-T E m) C t+BL mm (16) 
n 

and at finite B: HLee/ (cksP) — {i([jE{B — 0) + 
Bj M ,P E (B = 0) + BP M ]) - B 2 A MM } = B(-Tme) + 
B(-Tem) - (Tee) - B 2 {(T MM ) + hA MM /c}. 

The generalized kinetic energy Tem = —i[jET Pm] ap- 
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pearing in Eq. (|16f) is given by (in units hbc/K) 



at (in units hbc/K) 



M 



n 

^2 [ n Sn x (S n+ i x S n+2 ) + cycl. | /Vol 

n 

n 

-S„ x (S„+i x s„ +2 )| 

Vol 



where 'cycl.' means cyclic permutation of {n, n+1, n+2} 
- note that the ""-accent is permuted with the indices. 

(— Tem) changes sign under inversion of the S z - 
magnetization M and vanishes consequently for B = 0, 
as does the magnctothermal response in this case, com- 
pare (|16f) . The case of a finite B is therefore of interest. 
We will consider only the effect linear in B, since mag- 
netic fields amount normally only to small energies: 



8 



Leu ~ B——Lem 
oB 



B 



{-Tem(~[3M))-+Lmm 
(18) 

The EM-coefficient, therefore, consists of two parts. One 
is the -B-derivative of the zero-field Lem\ the other comes 
from the change of Je in magnetic field (namely, it ac- 
quires a term BJ) and is already linear in B - for small 
B. 



The quotient of Lme and L 



MM 



Q 



M 



me/Lmm 



Lme/Lmm 



(19) 



is called the (magnetic) thermopower. It has the advan- 
tage that geometric properties of samples cancel in ex- 
perimental studies and that the relaxation times cancel. 

To discuss Lme is superfluous, if one keeps the On- 
sager relation in mind. Anyway, if we follow the same 
steps as above, we arrive - with explicitly using OBC - 



Tme = -i[]M;PE 

—i 
Vol 



(n [s„ x S, i+ i, S„ • S„ + i 



(n + l)S„ x 



Sn+lj S„+i • S n+ 2 



71 S n+ i, S„ ■ S n+ i 



x S 



n+2 



= ^(n/2J x 2 (S„ - S n+X ) - (n+ 1)S„ x (S„ +1 x S„ +2 ) 
-nS n+2 x (S„ x S, i+ i)) 2 /Vol 

= ^n+l^n+2 ~ ^n+l^n+2 

n 

+ {n + 1) J z / Jx(Sn S n +l^n+2 + &n 

— nJz/Jx {S n S n+1 S^ +2 + 5^^1+1^1+2) 
+M/n B -S z Q •Vol/cl/Vol, 

which differs apparently from Tem- This does not 
contradict the Onsager relation because we find that 
(Tme) — (Tem) still holds. This may be seen by first 
using Eqs. (JSJ and ©: 



(P 



me 



EMI 



[[H,P m ].Pe] - [[H,P E ],P, 



M 



then adding a term which is zero by the cyclic property 
of the trace: 



= -([[1?,Pm],Pe] + [[P M ,P E ],H] + [[Pe.H].P m ] 



That this final expression is zero is just Jacobi's identity. 

Here it should be emphasized that (Tme) has non- 
negligible contributions from the boundary, and is - un- 
like (Tem) - sensitive to a change from PBC to OBC. 



VI. JORDAN- WIGNER TRANSFORM AND 
FREE FERMION MODEL 

For our spin Hamiltonian all quantities may easily be 
calculated when J z — via the Jordan- Wigner map- 
ping to a free fermion system. Under this condition 
the Hamiltonian is straightforwardly diagonalized by a 
Fourier transform with eigenvalues cos A: (setting J x = 1) 
following the Fermi-Dirac distribution (rife) = (c\c k ) = 

[l + exp(/Jcos k)] . 

The Drude weight entering (|18fl is given simply by 
-Lmm = (Tmm)c/H = j cosk(nk)dk/(2TTc)(n B c/% 2 ). 
At the same time Tem simplifies to (again in units 
fi B c/h): 

Tem = r^j J]] { — S^S^ +1 S n+2 — S n S^ +1 S^ +2 } 

n 

/dk . . 
— cos(2fe)n fc . 
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Following i|18[) we need to compute -right-hand side in 
units i± 2 B c/h - 



((-pM)(-T EM )) = ^§ E C0S ( 2A; ) ( n * ( n « - \ 



k , q 



{nkn q ) — (nk){n q ) holds for k ^ q, n 2 = n k and 
J2 k , q cos(2fc) (n k )(n q ) = 0. We find: 

((-/3M)(-Tem)) = -^]T[-cos(2fc)][<n fc )-K) 2 ] 



-0 



dk 



cos(2/c) 



27rc4cosh 2 [/3cos(fc)/2] 



(20) 



This result can be related with the temperature- 
derivative of the kinetic energy, 



(T BM {-f3M)) 



dT 



[T(-T MM }] , (21) 



where Tmm for the xy-Model is defined above. The re- 
lation 12111 is easily verified using a partial integration. 



VII. EINSTEIN RELATIONS 

Let us recapitulate Einstein's relation for diffusive 
transport in a metal: In a closed system the diffusive 
current jjj = — DVn, driven by a gradient in the particle 
density n, and the electrical current driven by an external 
potential, j = <jE, add to zero: jr> = —j- This condition 
yields Einstein's famous relation: 



a 

b 



Vn 
~~E 



(22) 



Here, in the case of a magnetic system, we are not in- 
terested in a variation of the electron particle density, but 
in a change in internal energy or magnetization. Hence, 
we replace n(x) — > A n where A n is either h n or /is 5^. 
Furthermore, the currents are not driven by a field E, 
but by a polarization Pp = ^ n nF n (times the lattice 
constant c), for a conserved F: [F, H] = 0. In our case 
the operator F is either the magnetization or the internal 
energy. 

We will also assume that there is diffusive transport for 
both the magnetization and the energy current. Thereby, 
we are led to the following definitions of the correspond- 
ing diffusion coefficients 

j M = DmmVM j e = Dee^E 

where M and E are the magnetization and internal en- 
ergy. 

Generalized Einstein relations can be derived using 
static response of VA\ n := (A, l+ i — A n )/c to a perturba- 
tion cPf, where c is the lattice constant along the chain: 



Vn/E 



\ X T (VA\ n cP F )/Vo\ 



where we perform a volume average. (This becomes 
mandatory because we consider the response to the cur- 
rent density rather than to the current at a fixed site.) 
Using the linearity of the isothermal susceptibility we find 
(by a discrete version of an integration by parts): 

c X T {VA\ n P F ) = J2{x T [A n+ i(mF m )]~ X T [A n (mF m )]} 

m 

= [( m + ^)x T {A n F m ) - mx T (A n F m ) 

in 

= J2x T (A n F m ) 

m 

and with the volume average on both sides 



^c X T (V^| n P F )/Vol = (3(AAAF)/Vo\ 



(23) 



We rewrite Q22t for the case of magnetothermal re- 
sponse, jM = Lme^T = Dmm^M = —jDi and find 



J M E 



L 



ME 



VAf 



VM 



D 



MM 



Dmm V/ 
- k B (3x T \VM\ n P E ) = 



k B T (VT/T) 

k B f3 2 



(24) 



Vol 



{AMAH) , 



for 



where we have assumed equal relaxation times r 
Lme — tLme an d the magnetization diffusion: Dmm = 
tDmm- (This assumption was only made for simplicity. 
At least, to our knowledge there is no reason why the re- 
laxation times should equal each other. However, the as- 
sumption is natural as in both cases the finite relaxation 
time comes from the fact that magnons lose momentum, 
and the physical processes responsible for that should be 
in both cases the same.) Eq. Q24|l leads to 



J ME 



D 



M M 



k B (3 2 
Vol 

Bk B p 2 
Vol 



(AMAH) a 
(AMAM) 



d k B p 2 
dB Vol 
Bk B p 3 



{AMAH) 



Vol 



(AM 2 AH){25) 



for small magnetic fields B, the structure of this equation 
is similar to Eq. (JTSJ. The first term of the right-hand- 
side of Eq. (HU is just Bx/T, where X = P{AM 2 )/Vo\ 
is the magnetic susceptibility. The second term of the 
right-hand-side of Eq. (|25|l simplifies to 



Bk B p 3 d 



{AMAM) = k B B/3 c 



d X 



Vol d/T ' " dp (3 

We therefore find the new relation: 
Lme 



B d 
TdT 



[T X ] 



D 



MM. 



= -X---[T X ] = -B d -* 
T TdT 1 AJ dT 



(26) 



Classically, the diffusion constant is D = v 2 t, where 
v is the velocity of the elementary excitations, here 
the magnon velocity. This leads at low-temperatures 
to a temperature-independent diffusion coefficient D = 
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D/t = v 2 and via (|26|l to vanishing magnetothermal re- 
sponse for T ~ * 0, whenever x(T) becomes constant for 
T -> 0. Setting A = M and F = M in the general 
Einstein relation yieldsii*i^ 



Lmm/Dmm 



X , 



(27) 



and with the choice A = H = F one may obtainS^i an 
analogous relation: Lee I Dee = cv- 

We will now introduce a formula which we believe to 
be approximately valid at small T. As stated above, 
the diffusion constant is known not to vary much near 
T = 0. Because of our restriction to low T we assume 
a temperature-independent diffusion coefficient Dmm 
(see discussion above). Then we can rewrite ()27(l as 
-§?L M m = -$f[x D MM] » D M M-§fX- Inserting this ex- 
pression into (|26|l we obtain 



ME) 



— B—Lmm- 



(28) 



In the case of free fermions this is just the result of Eq. 
(|2f \ . In the case of a finite interaction it seems to be 
correct in the limit T — > if we use data from Ref. 0. 
We provide tests of this relation in section ITXl 



VIII. AN EXACT EINSTEIN EQUATION 

In Sec lVIII we did set-up several versions of general- 
ized Einstein relations appropriate for magnetothermal 
response. Those relations may be viewed as a link be- 
tween the dynamical and static response theory, as they 
connect corresponding correlation functions by introduc- 
ing diffusion constants. 

Upon closer examination of the right-hand sides of the 
Einstein relations like namely, x, cy and dx/dT 

- we find that that these are static expectation values 
of products of AM- and AiJ-operators which could be 
generated by derivatives. Hence, it is an easy task to 
establish functional relations between the static correla- 
tions, e.g., 



dB 2 



d 

dT 



T' 



dT 



(TX) 



An intriguing question is, whether the analogous equa- 
tions obtained by switching between static and dynam- 
ical responses - x Lmm, etc. - could also be valid. 
For one of these relations, Eq. (|28|l . the range of validity 
is confined to the low temperature regime unless J z = 0. 
(A detailed discussion follows in Sec. lIXI ) One reason for 
the failure could be the fact that in this particular case 
one of the correlation functions does not reduce to a ther- 
modynamical expectation value. We therefore focus on 
a relation between the thermal and magnetothermal re- 
sponse coefficients - both mere thermal expectation val- 
ues. We claim that the relation 



dB 2 



-T EE ) 



dT 



[T 2 (T BM (-/?M)>] (29) 



is exact. The proof is provided in the appendix. This 
equation could be useful to derive an analytical solution 
for the magnetothermal response. 

The arguments of the proof rely mainly on the fact that 
one of the currents is a constant of motion. It is by the 
same line of arguments possible to show that Eq. Ij28(l is 
exact in the case where [H, Jm] = 0. Unfortunately, this 
is a very restricting constraint; only the xy- model and the 
Haldane-Shastry HamiltonianiS meet the requirement. 



IX. NUMERICAL RESULTS 

In the free fermion case Lem is easily accessible by a 
simple evaluation of the analytic expression 12U|) and a 
corresponding one for Lmm — ~(Tmm), using (|f 8fl . If 
J 2 / we have to compute a mere thermal expectation 
value in order to obtain Temi a task which is tractable 
to Monte-Carlo-simulations. Here, we use the stochastic 
series expansion (SSE), which is a Taylor expansion based 
Monte-Carlo- method^ 

We assume a Hamiltonian of the form H = 
J n h n , we may write the partition function Z = 



where the sum runs 



over all numbers m, all functions <p m : {1, . . . , m} — > N, 
and S^-eigenbasis-states a. The SSE-program now sam- 
ples the products of operators appearing in this expres- 
sion of Z with their relative weight factors. 

Normally, with MC-methods it is problematic to 
measure operators which are not diagonal in the S z - 
eigenbasis. But with SSE another class of operators is 
easily accessible: Following a standard procedure in sta- 
tistical mechanics, we have (h n ) = (dj n Z)/Z. For the 
SSE this means that we can measure any operator h n 
appearing in the Hamiltonian simply by counting how 
often it appears in the products sampled by our SSE- 
program (and dividing by J n ). So it is fairly easy to 
measure sums of products of parts of the Hamiltonian as 
S+S~ +1 , S~S^ +1 or S*Sn +v To make use of this fact 
in the given context, it is expedient to find an expres- 
sion for Tem which consists only of terms appearing in 
the Hamiltonian. Here we present one (Tem in units 
fiBc/ft): 



T, 



em 



^2{\jn,h n+ x] + [h n ,j n+1 }}/Vo\ 



n 

+JxJz ~ &n ^n+l^n+s] 

+J x Jz [S„Sn+n S^ +1 S n+2 — S n+1 S^ +2 \ I /Vol , 

which then allows to assess the magnetothermal coeffi- 
cient by the SSE. 

The Drude weight - which is the other input in Eq. 
(|18|l - may in principle be calculated by the Bethe-Ansatz 
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method. This was attempted by several authorSf2L2S but 
their calculations are still under discussion^ and so far 
there are no reliable results. We therefore use exact- 
diagonalization to obtain data for the Drude-peak. This 
has clearly the disadvantage that we cannot make any 
statements about the low-T and high-J z regimes where 
the convergence with Vol — ► oo is slow (cf. Ref. I24)) . 

In Figs. H and El we plot {{-(3M)T EM ) in units of 
J x [i B /h - a detailed discussion of units follows in the 
next section, compare Eqs. I|31|l and (|16|l - as a function 
of temperature for various interaction strengths. The 
data for J z ^ are obtained by MC-simulations and - 
in Fig. 0- by formula Eq. (gBJ and Eq. (fTHf> . 

The MC-results in Fig. [2 are clearly not good enough 
to determine exactly the T = 0- value. However, for 
Jz = Jx the data seems to extrapolate to 0.25, which 
is the Bethe-Ansatz result for the dimensionless T = 
Drude weight for the isotropic Heisenberg chain«i^ This 
result would imply via l|18|) a vanishing magnetothermal 
response for T — > and via l|28ll a vanishing T-derivative 
of the Drude weight, in contrast to previous results^ 

In Fig. [3 we present a comparison (at low T) between 
the MC-results for ((—(3M)Tem) and the results from 
equation (J2SJ, where we have used Bethe-Ansatz results 
(because of the finite size gap exact diagonalization pro- 
vides here no alternative) from Ref. |^] - which we prefer 
to Ref. [^because it is in agreement with Ref.|23|- for the 
temperature-dependent Drude weight Lmm- The agree- 
ment is very good for T < 0.2. As a consequence we may 
deduce from Eq. (gSJ that ({-(3M)T E m) = LmmH/c at 
T = 0, the linear magnetothermal effect as given by Eq. 
(|18|l vanishes consequently for T — > 0. 

The results for Lem and the thermopower (the prime 
denotes a derivative with respect to the magnetic field 
B) Q'm = L'me/Lmm are presented in Fig. 0] and Fig. [SJ 
respectively. Here we use exact-diagonalization, because 
for the computation of the Drude weight QMC-methods 
normally fail at higher T. (The standard proceeding is to 
extrapolate from the Matsubara frequencies to lo = as 
attempted in Ref. Q. This procedure becomes soon un- 
stable if T - and hence the spacing of the Matsubara fre- 
quencies - grows.) We use exact-diagonalization for the 
computation of (Tem) as well. Since (Tem) converges 
much faster with system size than the Drude weight, the 
finite size error is determined by the Drude weight, so us- 
ing QMC for (Tem) alone would not give better results. 
On the contrary, if we used MC-data we would introduce 
the statistical error. For the Drude peak we use Eq. 114|) 

with (u^)) = E En = E M^ k )( k \ J M^ En /z (z is 

again the partition function) for (Tem) we simply use 
the expression in Eq. (jT7Jl . The use of these expressions 
makes it possible for us to exploit translational symmetry 
which allows us to consider slightly larger systems than 
in Ref. |2J (namely, up to 18 sites). 

Because of the resulting finite size gap we can only 
discuss the high-T regime. (In the xy-model we find a 
vanishing thermopower at T = with a linear power- 
law.) 



For small J z we find a maximum at intermediate tem- 
peratures which decreases and is shifted to higher T as 
we increase J z . 

At higher J z the convergence becomes slower, so we 
cannot make any precise statements about the transport 
coefficients. However, in Fig.0]and Fig. [S]the curves de- 
crease with system size (for T sufficiently large). This 
behavior can be confirmed for all systems sizes (smaller 
than 18 sites). If we assume that the direction of con- 
vergence does not change for larger system sizes, our 
curves provide an upper bound to the exact results . But 
this would imply a reversion of the effect (negative ther- 
mopower) . While at J z — the effect is always positive - 
apart from a reversion due to the finite size gap, one finds 
Q'm < for a certain T-interval if J z is large enough. Be- 
cause of the slow convergence it is impossible for us to 
predict the precise location of this interval. Here one 
should emphasize that our conclusions were drawn from 
the inspection of relatively small systems, so we cannot 
really exclude that the observed reversion is only due to 
finite size effects. 

However, we see the reason for this reversion in the 
fact that the magnetothermal coefficient consists of two 
summands [see Eq. I|18|) ] which turn out to have opposite 
signs. These behave very differently under a change of J z . 
While the Drude weight dominates at J z = 0, it drops 
dramatically as J* is increased towards 1 (this was also 
reported in Ref. 22). The generalized kinetic energy does 
not decay so much, such that it becomes the dominating 
part. 



X. SIZE OF THE MAGNETOTHERMAL 
EFFECT 

Based on Eqs. and (|30l) and using the low B approx- 
imation in Eq. (|18f) we present in this section estimates 
of the size of the magnetic Seebeck effect. 



A. Dimensional analysis 

As we are considering effects linear in the magnetic 
field B we use / as a shorthand for J^. Using [Jm] = 
[/is]m/s and [Je] — Wm (compare Eq. (0 and (JBJ 
and Ijm/e] — [Jm/e]/™-) we niay decompose the re- 
sponse coefficients into dimension-full and dimension-less 
quantities 4 : 

T , _ k Bf i B fi B c 3 ~, 

L'em = — ^~^~h~^ x ^ 2 ^' EM ' 

where ~ denotes the dimensionless theory-result, c is 
the distance between neighbor sites, and hb is the Bohr 
magneton. Note that Lme has a factor 1/T with re- 
spect to Lem', hence a third factor J x (3 in (|30|l . In- 
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spection of that formula shows the units of Lme to be 
J(Tesla)- 2 m/(Ks 2 ). 

For the sake of completeness we state the analogous 
results 

T _ JxjlB MBC - 
LMM — 7 T—(JxP)^MM 

n h 

for the magnetic Drude weight and 



k B Jx JxC 



J EE 



(JxfifL 



EE 



for the thermal Drude weight. The units can be read 
off as [L M m] = J(Tesla)~ 2 m/s 2 and [Lee] = Jm/(Ks 2 ). 
Note that [L^] = [£y]s, see (j3J. 

We take the well studied one-dimensional Heisenberg 
chain S^CuOa as an example^ with the following pa- 
rameters c = 3.91 A, 25 J x /k B = 2.2 • 10 3 K A There are 
two chains in an area of 12.68 x 3.48A 2 iS This leads to a 
magnon velocity v s = J x c/h ~ 1.2 • 10 5 m/s. The mean- 
free path X S (T) extracted from a quasi-classical inter- 
pretation of thermal conductivity data& is strongly tem- 
perature dependent, it ranges from A S (200K) « 100 A 
to A S (50K) w 800 A. This results via A s = v s t in a 
relaxation time r(200K) ~ 0.8 • 10~ 13 s and r(50K) ~ 
6.6 • 10~ 13 s, in order of magnitude. 



D. Closed system 

Finally, we consider again the setting of a closed sys- 
tem. We are given a sample with open ends such that no 
current may flow. A temperature gradient should there- 
fore lead to a gradient in the magnetization. Using the 
Einstein relation for the magnetization current and Eq. 
IPS} we have 



VM = = 



D 



m .1/ 



\B = M is the magnetization caused by the presence of 
a magnetic field. Hence, dividing by M 

VM/M = -A[i n ^]VT 
or multiplying both sides with the system length 

AM/M = — -—. ^ Pixy] AT . 

Note, that AM/M is not directly a function of magne- 
tothcrmal coefficients. The magnitude of the dimension- 
less quantity ^[ln%] varies about 0.1. We therefore find 
with Jx/ks = 2200 K that the relative change in mag- 
netization - from one end of the sample to the other - 
should be a fraction of the order of 5 ■ 10 _5 K _1 x AT. 



B. Thermopower 

For the thermopower, we obtain 

n i _ L' ME _ k B , R \2^ME _ ^SA/ 

— 7 — ~r\ J xP) ~ = —j~WM ■ 

J-'MM Jx J->MM x 

Typically, k B /J x ~ lO^K" 1 (e.g., in Sr 2 Cu0 3 ) and 
at room temperature Q' M « 0.1. This yields Q' M w 
10- 4 K- X . 



C. Seebeck-effect 

The dimensionless _B-derivative of the magnetothermal 
response is (J x l3) 3 L' ME ~ 0.1, in order of magnitude, see 
Fig. H We therefore obtain from Eq. ij3U|) at T = 200 K 
for the magnetic particle current density a value of 

0.036 x B x VT moments/s 

per Tesla per (Kelvin/m), i.e., in Si-units. For a sample 
of length 1 cm = 10~ 2 m with a temperature difference of 
10 K and a field of 10 Tesla we have a current of 360 spin- 
1 moments (magnons) per second. A sample of Sr 2 CuC>3 
with a cross-section of 1 mm 2 contains 4.5 • 10 12 chains. 
It would induce a magnetic current of 1.6 • 10 15 moments 
per second. 



E. Peltier-effect 

We want to review now briefly the adjoint (Peltier) 
effect, described by the coefficient L EM . We now turn 
once again to our reference system Sr2Cu03. Here it 
is instructive to compare the energy currents driven by 
a temperature gradient and a gradient in the magnetic 
field. 

We set all quantities as above, with the exception of the 
values (J X P) 2 L' EM ~ 0.01, and {J X (3) 2 L' EE ~ 0.1, which 
we find more appropriate. We also do not compute the 
energy current density for one chain but per volume in 
Si-units. The result for the thermal response is 

je = 16 x VT J/(sm 2 ) 

and for the magnetic response 

3e = 3.3 • 10~ 4 x VB x B J/(s m 2 ). 

XI. CONCLUSION 

We have shown that non-trivial magnetothermal cur- 
rents may be induced in one-dimensional quantum-spin 
chains in an external magnetic field. We have argued 
that these effects might be especially important for the 
ID Heisenberg chain where the energy-current operator 
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commutes with the Hamiltonian and the magnetother- 
mal response diverges. We have presented estimates of 
the size for the current of magnetic moments induced by a 
thermal temperature difference for S^CuOa. We believe 
the size of the induced magnetic current to be sizeable 
and note that the observability of magnetic currents in 
magnetic insulators has been discussed recently^. 

The support of the German Science Foundation is ac- 
knowledged. 

XII. APPENDIX 
A. Proof of Eq. (1291 

1. Proof, first step 

We present the proof of Eq. (|29|1 in two steps. First we 
extend the Hamiltonian to comprise the static response 
to a thermal current: 

H{\,B)=H(B) + XJ E . 

We indicate the expectation value at nonzero A and B 
by indices. Moreover, we mean that these variables are 
zero, if we omit them in an index or an argument. 'A' 
means subtraction of the expectation value. 
We now show that 

{J M AH) X ,B = {JeAM) x , b . (32) 

We explicitly assume PBC - no polarization operators- 
and invoke the equation of continuity to obtain a local 
version of the above formula 

((Jm ~ Jm l )Ah m ) K B = -(i[H,AS z n ]Ah m ) x , B 
= {AS*i[H,Ah m ]}x,B = -(AS* n (J%- J%- 1 ))x,b. 

On the right-hand side of the above equation we perform 
a reflexion in space 

{(Jm ~ J^Ahm)^ = {ASl n (J E m ~ J E m+1 ))-^B . 
Translational invariance leads for 

a> n ,m '■— (JMAh m )\,B b n , m := — (J^AS^J-A.B 



to a„ jm = a„_ m and to 

Ofc — afc_i = bk+i — bk <=> afc — bk+i = a^-i — bk, 

where k — n — m. This last equation implies that 

c := a k - b k+ i 

does not depend on the index k. Our strategy is to show 
that c - a function on A and B - may be neglected with 
impunity in the thermodynamic limit. To this end we 
look at the following estimate: 

\c(B,X)\ = |&fc + i-afc| < min{|6 fc+ i| + \a k \}. 

k 

The fact that all correlation function of the form 
(AG1AG2) (as, e.g., a,k and bk) decay to zero when the 
spatial extension (i.e., the index k) goes to infinity leads 
to limvoWoo c = 0. Hence we find that in the thermody- 
namic limit afe = bk+\ => a k = X) therefore: 

(J M H)x, B /Vol = -(J e M)-x,b /Vol. 

The initial claim (|32|l follows - once again - by a reflexion 
in space on the second term. 



2. Proof, second step 

Having attained our first goal the further proceeding 
is standard. We apply consecutively derivatives with re- 
spect to A and B on our result and obtain: 

(JmJeAHAM) = -T^-(J E J E AM} B . 

U-D 

We note that the AH (resp., AM) might come from 
a derivative with respect to —d/d(3 = T 2 d/dT (resp., 
-Td/dB) and use (T EM (-0AM)) = -f3(JJ E (-/3AM)) 
and (—T EE ) = /3 2 (J E J E ). Our formula may then be 
rewritten as: 

^[T 2 {T EM (-(3M))] = T 2 ^(-T EE ) B , 

Q.E.D. 
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FIG. 1: Illustration of a quasi-one-dimensional magnetic in- 
sulator in the presence of an external magnetic field B, and a 
longitudinal temperature differential AT. A magnetic current 
jm, with the moments aligned along B is induced. 
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FIG. 2: Monte-Carlo-data (48 & 96 sites, PBC) for the di- 
mensionless ((—/3M)Tem) and various interaction strengths 
Jz/Jxx- The statistical MC-error-bars are given. 




FIG. 3: Monte-Carlo-data (48 & 96 sites) at small T for PBC 
in comparison (lines) with the estimates for {(—/3M)Tem) by 
Eq. Q280 and data from Ref. for the Drude weight. The 
cases J z = cos(7r/4) and J z = 1 are offset for clarity. 




FIG. 4: The linear part (in B) of Lem as a function of tem- 
perature, see Eq. 1180 . obtained by exact diagonalization for 
both ((-/3M)Tem) and the Drude weight. We plotted data 
for three system sizes (14, 16, 18 sites) with different line 
styles (long-dashed, dashed, solid, respectively). 




FIG. 5: The part of the thermopower 119H linear in B 
as a function of temperature, see Eq. 1181 with exact- 
diagonalization-data for ((— /3M)Tem) and the Drude weight. 
We plotted data for three system sizes (14, 16, 18 sites) with 
different line styles (long-dashed, dashed, solid, respectively) . 



